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1 Introduction 

Intermittency is one of the well-established routes to chaos, particularly in problems of 
fluid turbulence and in certain problems of reaction diffusion systemsfl, 2, 3, 4]. The 
appearance of alternate regions of regular and chaotic phases with the relative durations 
of two regions determined by the parameters of the system is the hallmark of this route. 
The saddle-node and the inverted pitchfork bifurcations [1, 2, 3, 5] form two important 
classes of bifurcations leading to intermittency. It is to be noted that while the saddle- 
node bifurcation is a typical codimension-1 bifurcation leading to type-1 intermittency, the 
inverted pitchfork bifurcation combines features of type-2 as well as type-3 intermittencies 
in the sense that an unstable branch of fixed points persists on the chaotic side of the 
bifurcation point. An exact renormalization group analysis of these bifurcations is also 
available in literature[3, 5]. 

Since quantum effects are ubiquitous, it is interesting to look into the effects of quan- 
tum fluctuations on intermittency. The introduction of quantum effects on a kicked rotor, 
that leads in the classical theory to the onset of chaos above a critical value of the param- 
eter, results in the suppression of chaos, similar to Anderson's localisation in the presence 
of disorder [6, 7]. The quantum fluctuations have been shown in our work to enhance chaos 
in both the types of intermittency (saddle-node and inverted pitchfork) referred to earlier. 
This enhancement is reflected in the three well-known features: (i) decrease of average 
staying time in the laminar (or regular) phase, (ii) increase of Lyapunov exponent, (iii) 
increased rate of diffusion in energy space. The phase space plots will also enlighten us 
about the quantum effects on the classical evolution. 

The squeezed coherent state formalism [8, 9, 10, 11, 12, 13, 14, 15] has been adopted in 
presenting a semi-quantal analysis of the theory. A similar analysis, followed also in other 
cases, has been seen to lead to enhancement as well as suppression of chaos depending on 
the examples chosen [9]. Starting with the Hamiltonian equations relevant to our problem, 
we have obtained certain model difference equations that give rise to intermittencies in 
the classical limit. This work is, we believe, the first attempt towards understanding the 
effect of quantum fluctuations on problems involving intermittency. 

Section II discusses briefly the theoretical background of the numerical calculations 
and Section III the numerical results. The last section concludes with the summary as 
well as the discussion of the results obtained. 



2 Formalism 

The squeezed state formalism[8, 9] stems from the time- dependent variational principle 

Ajdt(m\^--H\m) = o, (i) 

where ip(t) is chosen to be the squeezed coherent state represented by | i]){t)) = D(a)S( , y) \ 
0), D(a) = exp(<W — a*a), S(j) = exp[|(7at — 7*a 2 )] and | 0) represents the vacuum 



state. The pairs of canonical variables given by x 
Ax 2 = (x — (x)) 2 , Ap 2 = (p — (p)) 2 are obtained as 
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Clearly (x,p) represent the centroids of the wavepacket in the coordinate and the momen- 
tum space respectively while Ax 2 and Ap 2 give measures of the quantum fluctuations in 
the coordinate and the momentum space respectively. They satisfy the canonical equa- 
tions 
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where the Hamiltonian function H given by (ip \ H \ ip) can be expressed as 
H = ^ + V(x) + h(j- + 2n 2 <P) + {exp[^(^) 2 ] - l}V(x). 



(2) 



2.1 CASE I (Inverted pitchfork bifurcation) 

The mapping giving rise to the inverted pitchfork bifurcation in two dimensions corre- 
sponds to the periodically kicked Hamiltonian given by 

H = £ + ^ - |4 + ^) + *£<e + 3/tf) + 6 ^, E i(t _ „ T) , (3) 



T being the time interval between two successive kicks. 
Substituting for Ax 2 and Ap 2 from Eq. (2) we obtain 
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The Hamiltonian in Eq. (4) will yield the following difference equations after utilising 
Eq. (2) and going through algebraic manipulations: 

X n+ \ = X n + P n , (5) 

P n+l = P n + T[(eX n + (5Xl) + " ], (6) 



$ n+1 = $ n F n , F n = -^ + (n n + l) 2 (7) 

n 



and 



n n+1 = [IT n + n2 + -L]i + 2TA n+1 , (8) 

K F ™ 

where A n+1 = \[{e + 3/3X 2 ) + ™% X n = x n , P n = Pn T, $„ = 2 -^,U n = 2Tn n . 

The Eqs. (5) and (6) give, in the limit h — > 0, the classical inverted pitchfork bi- 
furcation. This is characterised by the fact that there exist three fixed points X = 
0, ±y( — f ), P = for e < 0, f3 > 0, the first one being stable, the other two being unsta- 
ble, going over to only one fixed point (unstable) at X = for e > 0, j3 > 0. The situation 
pertaining to e > 0, (5 < creates two stable fixed points at X = 0, ±w( — 4), P = and 
one unstable fixed point at X = 0, P = 0, while for e < 0, /? < the latter one becomes 
stable and the earlier ones vanish. The first case (/? > 0) corresponds to the subcritical 
and the second one (/? < 0) to supercritical bifurcations respectively in the classical case 
(i.e. h -> 0). 

2.2 CASE II (Saddle-node bifurcation) 

The saddle-node bifurcation corresponds to the Hamiltonian given by 

H = \( P 2 + Ap 2 ) - [{ex + ^- + ^f(2Px)} £ 8(t - nT). 
Substituting for Ax 2 and Ap 2 respectively we get 

H = [^ + |(_L + 47 r 2 0)] - [(ex + ^-) + *] £ tf(t - nT). (9) 

The corresponding difference equations are 

X n +i = X n + P n , (10) 

P n+1 = P n + T[e + /?X 2 + ^^], (11) 



$ n+1 = $ n F n , F n = — + (n n + l) 2 (12) 
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and 

n n+1 = [n n + n 2 + il +2TA n+1 , (13) 

V n r n 

where A n+ i = (3X n+ i (the symbols have their usual meanings). In the classical case 

(h —7- 0) these correspond to the existence of a stable fixed point X — — J(—jj),P = 
for e < 0,(3 > while for e > 0,(3 > no fixed point exists. This leads to chaos via 
intermittency. 

The phase space is divided into two regions, one bounded by X = ±X , P = ±P (-^o 
and Pq have suitably chosen magnitudes S and S' respectively) and called the laminar 
region, while the remaining region can be said to belong to the chaotic phase. Thus 
the journey in phase space will consist of regions of laminar behaviour interrupted by 
alternate regions of chaotic phase, characteristic of intermittency. The boundedness of 
the phase space ensures reentry of the phase points in the laminar region after a turbulent 
state and provided the number of iterations is large enough and the initial point suitably 
chosen, the reentry will be a uniform or white reentry. The average staying time in the 
laminar phase is calculated by assuming a uniform reentry which can also be simulated in 
the numerical calculation by starting with a uniformly distributed choice of initial points 
in the region. We have extended our calculations also to a non-uniform reentry. The 
numerical results remain qualitatively the same. 

3 Numerical Results and Discussions 

The signatures of quantum effects in classical intermittency through bifurcations, albeit 
in the semi-quantal domain, were investigated by (i) considering the phase space plots and 
(ii) probing into other hallmarks of chaos (given later) relevant to two maps exhibiting 
saddle-node and inverted pitchfork bifurcations respectively. The exit of the phase space 
points starting from the neighbourhood of stable fixed point becomes, as is evident from 
the Figs. 1-3, progressively more and more rapid with increase of Ti in both the maps, 
thus supporting our conjectures about the quantum enhancement of chaos. The quantum 
tunneling effects are demonstrated by the plots (Fig. 2), where the system is a single 
particle periodically kicked by a potential which is a double well with a single hump in 
between, with the initial point chosen to be close to any one of the stable minima. The 
other plots given subsequently (Figs. 4-11) corroborate our view in this regard. 

The relevant features of chaos, namely (i) the average staying time in the laminar 
phase, (ii) the Lyapunov exponent, (iii) the energy diffusion were determined for both 
the above maps. Since we have endeavoured to probe into the situation on either side 
of the bifurcation point, the magnitudes of the parameters e, (3 were chosen accordingly. 
For the inverted pitchfork bifurcation the parameters chosen were (i) e > 0,(3 > 0, (ii) 
e < 0, (3 > and (iii) e > 0, (3 < respectively. The first case corresponds to the existence 



of an unstable fixed point at X = 0, P = 0, while the case (ii) corresponds to the presence 
of three fixed points of which the pair at X = ±v/( — 4), P — are unstable and the one 
at X = 0,P = is stable (clearly here e = is the local bifurcation parameter). The 
case (iii) corresponds to supercritical bifurcation and we have included this case as well 
to show the effect of two stable fixed points at X = ±v/( — -§), P = on the relevant 
parameters both in the presence and the absence of quantum effects. 

The numerical calculations have been performed in the following way: (a) In the 
calculation of the average staying time in the laminar phase an ensemble of initial points 
(assuming both uniform as well as non-uniform distribution) surrounding the stable fixed 
point was chosen and for each of these the number of iterations required for exit out of the 
laminar region (S = S' = 0.1 here) was calculated. The average was calculated for each 
value of ft, (b) for calculation of the average Lyapunov exponent A, a pair of closely lying 
initial points (with separation, say, do) was chosen near the stable fixed point. After every 
two iterations the distance d between the two points was calculated and the expression 
Ai = | In j- determined numerically. The distance was normalized to the initial value and 
the process repeated till the evolved points left the laminar region. The resultant average 
Lyapunov exponent A over the entire laminar region was then computed, (c) lastly, for 
the estimate of the energy diffusion, average value of p 2 was calculated for a choice of the 
initial points lying close to the stable fixed point for each value of %. In all these cases, 
the calculations were repeated for different choices of the parameters and initial points 
respectively. Since the results obtained were similar, only a typical set has been displayed 
in the manuscript. 

The influence of the stable fixed points is evidently manifested, at least for low quan- 
tum effects by (i) the increase of average staying time in the laminar phase, (ii) the 
decrease of Lyapunov exponent, (iii) the resistance to the energy diffusion compared to 
that in the absence of the fixed point. For example, the parameters (a) e < 0, (3 > and 
(b) e > 0, j3 < correspond to the existence, while the parameters (c) e > 0, /? > pertain 
to the absence of the stable fixed points in the inverted pitchfork bifurcations. For the 
saddle-node bifurcation we have investigated the case e < 0,/3 > as well as the case 
e > 0, j3 > 0. The first one corresponds to a stable fixed point at X = — J( — |), P — 
while there is no fixed point on the other side of the bifurcation point e > 0,/3 > in 
the second case (The case e > 0, (3 < is uninteresting in this case). A range of possible 
values of e, j3 was covered but since the results were found out to be qualitatively similar, 
the results for a typical choice of parameters have been displayed here. 

It should be noted that although the effects of quantum fluctuations on classical dy- 
namics have been investigated in some standard Hamiltonian maps [9] no such attempt 
has been made earlier in problems involving chaos via routes of intermittency through 
bifurcations. The semi-quantal formalism gives rise to a set of four coupled equations 
involving x, p, (Ax) 2 , (Ap) 2 . The effect of the quantum fluctuations (Ax) 2 , (Ap) 2 on the 
evolution of the wavepacket leads for a non-compact phase space relevant to our problem 
to an enhancement of chaos. Analogous conclusions have been reached in some earlier 
works[8, 9, 10]. 



In both the maps we have dealt with the situations corresponding to the presence of 
unstable fixed points (or no fixed point in the case of saddle-node bifurcations) as well 
as stable fixed points lying on the opposite sides of the bifurcation points respectively. 
The variations of the relevant parameters with the increase of quantum effects are shown 
graphically. The numerical calculations yield the following features corresponding to both 
the maps. 

If we start with initial points lying close to stable fixed points (in the basin of attraction 
of the fixed point) (which are (0,0) for e < 0, /3 > 0, and (±\J~( - -§), 0) for e > 0, /3 < 

in the case of the first map and (— J( — §),0) for the second map respectively), the 
average staying times in the laminar phase will be, for low h, greater than those in the 
case corresponding to the absence of any stable fixed point (e > 0, j3 > or e < 0, (3 < 
in both the maps). This is quite expected as a consequence of the presence of the stable 
fixed point in an almost classical domain. Contrary to classical dynamics, where the 
initial points close to the stable fixed point remain always in its neighbourhood, there is 
a wandering of the points away from the stable fixed point for large quantum effects, i.e., 
for large values of Ti. This is also amply manifested in the phase space plots shown (Figs. 

1-3). 

The classical effects are seen to be offset by the introduction of quantum fluctuations 
which decrease the average staying times in the laminar phase leading to practically 
identical values for large values of % (Figs. 4 and 5). In respect of the average Lyapunov 
exponent A, the observations are similar except that A increases with the increase of ft 
(Figs. 6 and 7) and for low h it is smaller in the presence of a stable fixed point. Regarding 
the energy diffusion, the presence of the stable fixed point decreases the rate of diffusion 
compared to the case of no-stable fixed point for low value of h. Moreover the energy 
remains bounded within certain limits leading to a certain kind of localisation in the case 
of the initial points chosen close to the stable fixed points at least for low quantum effects 
(Figs. 8 and 10). For higher quantum effects the energy diffusion becomes much more 
rapid (Figs. 9 and 11). However, for large quantum effect, contrary to the classical case, 
no simple scaling law [5] regarding the variation with e has been seen to exist. Though 
we have confined ourselves to the laminar region, the situation is expected to be similar 
for large % in the case of the turbulent region which lies outside the gate bounding the 
laminar phase. A different choice of initial X and P, for example, both lying in the region 
of classically unstable fixed points has been seen to display a similar behaviour. 

The enhancement of stochasticity (reflected in the decrease of the average staying 
time in the laminar phase, increase of the average Lyapunov exponent as well as the rate 
of energy diffusion with the onset of quantum effects) can be qualitatively explained by 
condsidering predominantly the non-compactness of the phase space and also partly by 
taking into account the broadening of the wavepacket as well as the quantum mechanical 
tunneling. The significance of the typical quantum mechanical features is revealed in 
the graphical plots where the preponderance of the quantum mechanical effects is clearly 
visible even when the initial points are chosen to lie in the neighbourhood of the clas- 
sical stable fixed points. Similar quantum mechanical enhancement of chaos has also 



been observed if the initial points are chosen to lie near the classically unstable fixed 
points (although not shown in the plots displayed here). It is, therefore, clear from our 
results utilizing the squeezed coherent state formalism that the quantum tunneling and 
wavepacket spreading effects combined with the non-compactness of the underlying phase 
space dominate over the chaos-suppressing interference effects [16] for the type of bifur- 
cations investigated. Similar results have also been obtained by the earlier workers[8, 9], 
albeit in other contexts. 

The squeezed coherent state formalism thus enlightens us about the effects of quantum 
fluctuations on the problems of intermittency involving bifurcations of the saddle node 
and the inverted pitchfork types respectively. This semiquantal formalism admits of 
applications to other problems of classical chaos. 
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Figure Captions 

Fig-1. Phase space plots for the inverted pitchfork bifurcation (e = — 10~ 4 , (3 = 1.0, 
(f) = 0.5 (initial value), n = 0.0 (initial value) and T = 1.0) with the initial point close to 
the stable fixed point X = 0, Po = for fi = (left figure), = 10~ 6 (middle figure), = 1.0 
(right figure). The enhancement of chaos with the increase of quantum effects is reflected 
in the progressively rapid exit of the phase space points with the increase of h. 

Fig-2. Phase space plots for the inverted pitchfork bifurcation (e = 10 -4 , (3 = —1.0, 
<fi = 0.5 (initial value), n = 0.0 (initial value) and T — 1.0) with the initial point close to 
the stable fixed point X — ±0.01, P — for h — (left figure), = 10~ 6 (middle figure), 
= 1.0 (right figure). The enhancement of chaos with the increase of quantum effects is 
reflected in the progressively rapid exit of the phase space points with the increase of fi. 

Fig-3. Phase space plots for the saddle node bifurcation (e = — 10~ 4 , (3 = 1.0, 4> = 0.5 
(initial value), n = 0.0 (initial value) and T = 1.0) with the initial point close to the 
stable fixed point X = — 0.01, Po = for fi = (left figure), = 10~ 6 (middle figure), 
= 1.0 (right figure). The enhancement of chaos with the increase of quantum effects is 
reflected in the progressively rapid exit of the phase space points with the increase of %. 

Fig-4. Plot of average staying time against log 10 (h) for the maps exhibiting inverted 
pitchfork bifurcation. The three graphs correspond to the cases a) e > 0,/3 > 0, b) 
e < 0, (3 > and c) e > 0,(3 < respectively. The parameters chosen are e = 10 -4 , 
(3 = 1.0, (f> = 1.5 (initial value), n = 0.0 (initial value) and T = 1.0 respectively. The 
laminar region is given by Xq — Po = 0.1. The phase space is chosen as a square of 
size 1.0. The decrease of average staying time with the increase of quantum effect is 
reflected in the figure. The results shown here are for a uniform distribution of initial 
points although they were seen to be similar when this distribution is non-uniform. 

Fig-5. Plot of average staying time against log w (h) for the maps exhibiting saddle- 
node bifurcation. The two graphs correspond to the cases a) e > 0,(3 > and b) e < 
0,f3 > respectively. The parameters are identical with those in the first case. The 
laminar region and the phase space are chosen as before. The results shown here are for 
a uniform distribution of initial points although they were seen to be similar when this 
distribution is non-uniform. 

Fig-6. Plot of log 10 (average Lyapunov exponent) against log 10 (/i) for the maps 
exhibiting inverted pitchfork bifurcation. The three graphs correspond to the cases a) 
e > 0, (3 > 0, b) e < 0, (3 > and c) e > 0,(3 < respectively. The initial points are 
chosen to lie in the neighbourhood of the stable fixed point for the last two cases. The 
parameters chosen are the same as in the earlier cases except that the initial value of (f> is 
0.5. do (see text) is chosen close to 10~ 9 . 

Fig-7. Plot of log 10 (average Lyapunov exponent) against log 10 (h) for the maps ex- 
hibiting saddle-node bifurcation. The two graphs correspond to the cases a) e > 0, (3 > 0, 
b) e < 0, (3 > respectively. The initial points are chosen to lie in the neighbourhood of 
the stable fixed point for the second case. The parameters chosen are the same as in Fig. 
6. 
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Fig-8. Plot of log 10 (energy) against the iteration number a for the maps exhibiting 
inverted pitchfork bifurcation in the cases a)e>0,/3>0, b)e<0,/3>0 and c) 
e > 0,/3 < respectively. These plots refer to low quantum effect (h = 10~ 10 ). The 
values of the parameters are chosen as in Fig. 7. The initial points are chosen to lie in 
the neighbourhood of stable fixed points for the last two cases. 

Fig-9. Plot of log 10 (energy) against the iteration number a for the maps exhibiting 
inverted pitchfork bifurcation in the cases a) e > 0, j3 > (the solid line), b) e < 0, [3 > 
(the long-dashed line) and c) e > 0, j3 < (the short-dashed line) respectively. These 
plots refer to large quantum effect {Ti = 1CP 2 ). The values of the parameters are chosen 
as in Fig. 7. The initial points are chosen to lie in the neighbourhood of stable fixed 
points for the last two cases. The first two cases (the solid line and the long-dashed line) 
give almost identical results for this value of U and hence are indistinguishable in this 
plot. 

Fig-10. Plot of log 10 (energy) against the iteration number a for the maps exhibiting 
saddle-node bifurcation in the cases a) e > 0,/3 > and b) e < 0,/3 > respectively. 
These plots refer to low quantum effect (h = 10~ 10 ). The values of the parameters are 
chosen as in Fig. 7. The initial points are chosen to lie in the neighbourhood of the stable 
fixed point for the second case. 

Fig-11. Plot of log 10 (energy) against the iteration number a for the maps exhibiting 
saddle-node bifurcation in the cases a) e > 0,/3 > and b) e < 0,(3 > respectively. 
These plots refer to large quantum effect (h = 10~ 3 ). The values of the parameters are 
chosen as in Fig. 7. The initial points are chosen to lie in the neighbourhood of the stable 
fixed point for the second case. 
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